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We investigate possible nontrivial phases of a two-subband quantum wire. It is found that inter- 
and intra-subband interactions may drive the electron system of the wire into a gapped state. If 
the nominal electron densities in the two subbands are sufficiently close to each other, then the 
leading instability is the inter-subband charge-density wave (CDW). For large density imbalance, 
the interaction in the inter-subband Cooper channel may lead to a superconducting instability. The 
total charge-density mode, responsible for the conductance of an ideal wire, always remains gapless, 
which enforces the two-terminal conductance to be at the universal value of 2e^//i per occupied 
subband. On the contrary, the tunneling density of states (DOS) in the bulk of the wire acquires a 
hard gap, above which the DOS has a non-universal singularity. This singularity is weaker than the 
square-root divergency characteristic for non-interacting quasiparticles near a gap edge due to the 
"dressing" of massive modes by a gapless total charge density mode. The DOS for tunneling into 
the end of a wire in a CDW-gapped state preserves the power-law behavior due to the frustration 
the edge introduces into the CDW order. This work is related to the vast literature on coupled ID 
systems, and most of all, on two-leg Hubbard ladders. Whenever possible, we give derivations of 
the important results by other authors, adopted for the context of our study. 
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I. INTRODUCTION 



From. a theorist's point of view, electrons in quantum wires should provide a simplest realization of a Luttinger 
liquidEJH. Indeed, as the motion is confined in the direction transverse to the axis of a wire, the system is effectively 
one-dimensional; the electron-electron interaction is strong enough (typically, of the order of the Fermi energy) for the 
interaction effects not be washed by the temperature; and the state-of-the-art wires (at least the semiconductor version 
of them) are clean enough for disorder effects to be sufficiently weak. A crucial experimental test for the existence of 
the Luttinger-liquid state in quantum wires would be provided by tunneling into a wire, either in the middle or into 
the end. As is well-known, the tunneling density of states (DOS) of a Luttinger liquid reveals a pseudogap behavior, 
i.e, it is suppressed at energies close to the Fermi-energy, which should result in a power-law bias-dependence of the 
tunneling current, and in a power-law temperature dependence of the Ohmic conductance. Finite-bias and finite- 
temperature measurements have already been performed on a very special realization of quantum wires - carbon 
nanotubes, and observed non-linear current-voltage dependences were interpreted in terms of the Luttinger-liquid 
theoryHu. Features of resonant tunneling, characteristic for a Luttinger liquid, have recently been observed on GaAs 
quantum wireaJ prepared by cleaved edge overgrowth technique. Also, a LuttiiJger-liquid behavior has been reported 
in tunneling into InSb wires naturally grown in a porous material (asbestos)EL Tunneling pseudogap oLa quantum 
wire has been described in terms of a Luttinger-liquid model both for a single-l3 and multi-subband wirds, the latter 
system exhibiting a smooth healing of the pseudogap as the number of the occupied channels increases. In anticipation 
of more and better controlled tunneling experiments on quantum wires, and also from a general point of view, we 
would like to ask if there are any processes which could open a true gap, rather than a pseudogap, in the electron 
spectrum of a wire, and if yes, what are the properties of the corresponding gapped phases. 

To this end, we consider in this^paper a two-subband quantum wire, having in mind semiconductor nanostructures 
studied recently in, e.g., Refs.u" EiJ. To some approximation, this system is similar to other two iHiportant classes 
of ID two-band systerns-siLidied extensively over the last few years, i.e., two-leg Hubbard ladderstJiiS and (single- 
wall) carbon nanotubeaifTo (for an account of earlier results on coupled ID systems, see, e.g.pj). Studies of Hubbard 
ladders identified inter-subband scattering processes capable of driving the system into a gajjaed state. Phase diagrams 
of a generic ladder, containing a multitude of gapped states, were constructed in Refs.tjtj. A similarly formulated 
problem, with appiications to a ID system with electron and hole bands (valence-fluctuation problem), was investigated 
some time ago inc3. 

The goal of the present paper is two-fold. First of all, we would like to understand which of the gapped phases, 
found in Hubbard ladders, have a chance to occur in quantum wires. The main difference between these two systems 
is that the Coulomb interaction in wires is (i) (supposed to be) purely repulsive; (ii) relatively long-ranged (even 
in the presence of a metallic gate); (iii) relatively well-known at distances larger than the lattice spacing (which is 
the range relevant for quantum wires) ; this imposes constraints on the choice of coupling constants for the Hubbard 
model. Also, because the electron wavelength in wires is larger than the lattice spacing, Umklapp scattering is 
unimportant. All these constraints reduce the variety of possible gapped states to (a) inter-subband charge-density 
wave (CDW), and (b) superconducting state. [We will come back to a more detailed description of these states shortly] 
The main difference between carbon nanotubes and quantum wires is that the former, because of its special crystal 
structure, has two conducting subbands with commensurate Fermi-momenta. Although, as we will show, a quantum 
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wire with nominally different subband Fermi-momenta may be driven into the commensurate state by inter-subband 
backscattering, this process occurs in a competition with other processes and requires special analysis. 

Having determined which gapped phase can in principle occur in a quantum wire, we focus on the calculation of 
measurable quantities in each of these phases, which is the second goal of the present paper. Our main emphasis is 
on the tunneling density of states, which turns out to exhibit an unusual threshold behavior due to coexistence of 
gapped and gapless modes and also be sensitive to the presence of open boundaries. In addition, we consider the 
two-terminal conductance both in the absence and in the presence of impurities. 

Although this paper is not stipposed to be a review, we present, when possible, derivations of important results by 
other authors, e.g., Refs.c3t3ca, adopted to the context of our study. Hopefully, this would help a reader, who is not 
an expert in the field, to understand connections between different approaches. 

Having formulated the goals and scope of this paper, we now return to a generic two-subband quantum wire 
with incommensurate Fermi-momenta kip and k2F in subbands 1 and 2, respectively. In the basis of occupied 
transverse states, this becomes the problem of two Luttinger liquids coupled by inter-subband interactions. To 
understand possible phases of such a system, one should consider all possible scattering processes involving electrons 
from different subbands, i.e., forward scattering, backscattering, and "Cooper scattering" (cf. Figures 0,111). Forward 
scattering simply renormalizes parameters of Luttinger liquids formed by electrons of each subbands but does not 
result in new phases, although it does change the conditions for occurrence of new phases. (To be more precise, forward 
scattering between Luttinger liquids in different subbands is responsible for the crossover into the Fermi-liquid state, 
but this crossover occurs smoothly as the number of channels increases) . 

The momentum transfer in a backscattering event involving electrons of different subbands is equal to kip ± k2F- If 
\kip — k2F\ ^ T/imiii{vip^ V2f}^ then there are no final states available for electrons involved in such a process (here 
Vip and V2F are the Fermi velocities in the two subbands). Thus, if the temperature T is low enough, interchannel 
backscattering is forbidden. However, it may become energetically favorable for a system to equalize the charge 
densities, and hence the Fermi momenta, of different subbands. This may occur if the Fermi-momenta difference is 
small enough and the amplitude of backscattering is large enough. After the densities are adjusted, backscattering 
becomes possible. As a result, inter-subband charge-density wave (CDW) phase may be formed, in which charge 
densities of the subbands form a staggered pattern, see Fig. a Similarly to classical charge-density waves, this phase 
is very sensitive to a random potential, resulting in pinning of the CDW and strong suppression of conductance with 
disorder. 

The "Cooper" scattering event, on the other hand, always conserves momentum and energy. In this process, two 
electrons starting in, e.g., subband 1 with momenta kip and —kip, scatter on each other and end up in the other 
subband, also with opposite momenta k2F and — fc2_F- "Cooper scattering" can be considered as formation of a 
fluctuational Cooper pair in one of the subband followed by its tunneling into the other one. When kinetic energy 
gain due to such tunneling overcomes the Coulomb repulsion, the wire is in a Cooper (or superconducting) phase. 
This phase is characterized by locking of fluctuating charge currents in different subbands to each other as well as by 
spin gaps in each of the subbands. The Cooper phase is favored when Fermi-momenta imbalance is largest, i.e. when 
the second subband just starts to fill up. Disorder has a less pronounced effect on Cooper phase than on CDW one, 
similarly to what happens in higher dimensions. 

It is important to emphasize here that inter-subband backscattering and Cooper scattering block only modes 
corresponding to relative charge- and spin-excitations, but leave the center-of-mass charge mode free. As a result, the 
conductance remains at the universal value of 2e^/h per occupied subband in clean CDW and Cooper phases. 

Despite being ideal conductors, both the CDW and Cooper phases are characterized by the truly gapped behavior of 
the tunneling density of states at energies below corresponding gaps. This is so because a ID electron is a convolution 
of charge and spin collective excitations, and if some of these excitations acquire a gap, the entire electron acquires it 
as well. 

Somewhat surprisingly, we find that tunneling into the end of the CDW wire is quite different. A tunnel barrier at 
the end of the wire distorts charge-density wave profile and creates a static semi-soliton. This allows tunneling into 
the end to oc cur e ven at energies below the bulk CDW gap (in the lowest order in the barrier's transparency). 

In Section VI B we consider, for illustrative purposes, the density of states-pLlhe "Mott phase" , which occurs 
in a single- or multi-subband wire subject to an external periodic potcntiaEjC In the case of a wire, formed 
in a semiconductor heterostructure, this potential may be provided by an additional electrostatic gate of periodic 
shapeEil. Varying the potential applied to this gate, one can tune electrons of the wire into the half-filling condition 
(one electron per unit cell of the periodic potential). Unlike the two strong-coupling phases mentioned above, the 
Mott phase, which is described by the half-filled Hubbard model, does not conduct current because its total charge 
fluctuations are gapped by the external potential. 
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II. HAMILTONIAN OF A TWO-SUBBAND QUANTUM WIRE 



A. Classification of scattering processes 

Electrons in a quantum wire are described by the following Hamiltonian 

+ IJ2J dVdV'C/(r-r')*Ur)*l,(r')*s'(r')*s(r), (1) 

where s is the spin index, 14on/ (r_L) is the confining potential in the transverse direction, and U(r) is the electron- 
electron interaction potential. The Fermi-wavelength of electrons is assumed to be much larger than the lattice 
spacing of the underlying crystal structure. Because of t hat, w e do not consider umklapp processes, in which electron 
momentum is transferred to the lattice (except for in Ch.VIB). Hamiltonian (|l|) is Galilean-invariant, and hence our 
subsequent calculations have to preserve this invariance as well. We will return to this important point later on in 
our discussion. 

If the chemical potential in the leads is such that only two lowest subbands of transverse quantization are occupied, 
the electron wavefunction is given by 

2 

^-,(0 -^0„(ri)V'„.(x), (2) 

n=l 

where (/)„(rj^) are the orthogonal wavefunctions of transverse quantization, chosen to be real. In this basis, the kinetic 
part of Hamiltonian (mj becomes 



(x) (3) 



where e„ is the energy of the n-th transverse subband. 

To describe low-energy excitations in the n-th subband, we expand the longitudinal part of the ^'-operator, tpsi^x), 
in terms of right- and left-moving excitations, residing around the Fermi-points of the n-th channel: 

i^Ux) = i?„,(a;)e^'="-" + L„,(x)e-*^-"-^ (4) 

In this representation, the interaction (four-fermion) part of Hamiltonian (|l|) reduces to a sum of two terms. The 
first one, Umtra, describes the interaction of electrons within the same subband, and contains forward and backward 
scattering processes. The second one, Umter, describes the inter-subband interaction. It splits naturally into forward 
(U^), backward (U^), and Cooper (U^) parts. 

U^nter = + + . (5) 

The momentum transfer between subbands is limited by the width of the Fermi distribution function, and is therefor 
small at low temperatures. Forward scattering involves no momentum transfer between subbands (cf. Figj^). This 
process is also an example of a direct process, in a sense that electrons stay in the same subband, as is evident from 
the explicit expression for 



X [RlA^')Rr,^Ax') + LlA^')L„^Ax% (6) 

where the direct matrix element is given by 

m\-^-\x~x') = I [/(r-r')0,^,(r^)^^(rl); (7) 

and ^, = J dz J dz' . 
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By "backward scattering" , we understand processes with a non-zero momentum transfer Sk = kip ± k2F between 
subbands. These processes can be divided further into direct and exchange parts 

[/^ = C/f + t/f. (8) 

In an exchange process [see Figs.^,||, electrons change subbands. Two parts of backscattering can be written as (for 
the sake of brevity, we omit here the x — x'-dependence of the matrix elements) : 



n^TTi s,s' 

{R L)e-2»(fc™F-fc„F)(:i^+x')] (9) 



and 



= E E / Mi-^mRlix)Rr.Ax')RlA^')Rrasia 

n^m s,s' 

+ {R ^ L)e-*('='"^-'="'=')(^-^')} 

+{i?L(^)^n.'(:j')4.s'(^')ims(a;)e'("'"-+'="-)("^"') + {R ^ L)e~'^ik^-+k^-)i—')}l (10) 
where the exchange matrix element is 

Mt'^Hx-x')^ f C/(|r-r'|)0„(r^)0„(rl)0„(rl)0„(r^). (11) 




FIG. 1. Example of inter-subband forward scattering. Filled (empty) circles denote initial (final) states of electrons. Dashed 
lines with arrows indicate "direction" of the scattering. 




FIG. 2. Example of inter-subband backscattering: (a) direct, (b) exchange. Notations as in Fig.|l| 



FIG. 3. Examples of inter-subband Cooper scattering. Notations as in Fig.|l| 

Momentum conservation requires that the energy of at least one of the states, involved into direct backscattering, 
should be far away from the Fermi energy, which forbids this process at not too high temperatures (T <ti \kmF — 
knpl min {vnF,VmF})- This is reflected in the presence of the exponential factors in front of the fermion operators in 
Eq.(H), which oscillate rapidly as functions of (x + x'). This restriction can be lifted though, if the system prefers to 
gain energy from backscattering by equalizing the subband densities, so that kmP — knP- 

Finally, we call "Cooper scattering" (Figj^) a process in which two electrons with zero total momentum (a fluctua- 
tional Cooper pair) hop from, e.g., channel m , into channel n, so that the total mcmpritum Q = —km + fcm — rAr=f 
— fc„ + fcn = is conserved. This process is also referred to as "Josephson couphng"E§E3, or as "gooTTTr process"ll§ll3. 
The Hamiltonian of Cooper scattering is given by 

By construction. Cooper scattering is of the exchange type. 

In what follows, we use the following abbreviations: forward scattering = FS, direct backward scattering = dBS, 
exchange backward scattering = xBS, Cooper scattering = CS. 

For a generic situation of fc„ ^ k^, the only momentum-conserving inter-subband scattering processes are FS, xBS, 
and CS. The amplitudes of these processes depend on the ratio a/d , where a is a typical transverse size of the wire 
(which determines the spatial extent of 0„(rj^)) and d is the interaction range. In the limit of a/d — > 0, the interaction 
potential can be taken out of integrals (0) and (|ll|), upon which remains finite, whereas vanishes. It can 
be readily shown that for finite but small ratio a/d, the exchange matrix element is small: Mx ~ {a/d)'^Md. The 
long-range interaction thus discriminates against exchange processes. If (as it is most often the case) a wire is formed 
by means of a gate deposited over the 2D heterostructure, d is given by the distance to the gate, which screens the 
Coulomb interaction in the wire (see Fig. ^). Typically, a/d = 0.1 — 1. 




FIG. 4. Schematic view of a gated wire. The wire of typical transverse size a is separated by distance d from the metallic 
gate. Distance d determines the range of interaction among electrons inside the wire. 



B. Bosonized form of the Hamiltonian 



We use the conventional bosonization procedure, in which 

= ^=e'V^('^"=-^"=), (13) 

V 27ra 

LU^) = ^=e-'^^('^"=+^"=), (14) 

V 27ra 

and short-range cut-off a ^ kp^. Boson fields (fns and 6ns with n = 1,2; s = ±1, are decomposed into charge- (p) 
and spin- (cr) collective modes 

0ns^ -^{Onp + S0na)- (15) 
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Parts of the Hamiltonian, corresponding to free motion and intra-subband interactions {Hq + Uintra), take the well- 
known Luttinger-liquid form: 

Hnp = 77 / dx{VnpKnp[dxOnpY + -fT^ {dxfnpf} , (16) 
^ J J^np 

If V 

+lSy^ / (17) 

which describes independent charge- and spin-density excitations {Hnp and Hna-i respectively). The cosine term in 
Eq. ( |l7|) is due to backscattering within a single subband. Explicit expressions for the Luttinger-liquid parameters will 
be discussed later. 

Upon bosonization, the three types of intersubband interactions take the following form : 



^ — I dxdx(pipdx(p2p; (18) 



^ J dxdxipipda:(p2p; 



Uf = „ ^ / dxcos 



C/f / dx 



^{tpip - ip2p) + 2{kiF - fc2_F)a 
X C0s[V2TTipi„] COs[/2Tiip2a]] (19) 

/" , / ^1 + ^2 (^Q^ip^^g^ip^^ ^ da:ipiadx^2a) + — — {dxdlpdx02p + da:0i„dx02a) 

V 77 77 

2 

C0s[\/277(6'ic, - 02a)]{ibl + ^2) COs[\/277l^io.] COs[\/277l^2cr] - 



77^ 



(61 - 62) sin[V2^(pi<,] sin[V2^^2^]} ) ; (20) 



t/ = ^ , , / da;{isp COs[V277(0lp - 6'2p)] COs[V277(y5lCT] COs[V277(y52CT] + 

277^a^ / 



-Kitp cos[%/2^(6'ip - 6l2p)] ( cos[\/2^(6li„ - 02^)] - sin[\/2^^i^] sin[\/2^(y92„]) }. (21) 
The corresponding amplitudes are given by 

/o^y dxM^/^\x), 

fbs = J dxAI^^^x) cos[{kiF + k2F)x], 

6i,2 = J dxAfJ^^^(a;) cos[(fcii^ T ^2;^)^;], 

tgp ~ J dxAI^^'^^{x) cos{kipx) cos{k2Fx), 

ttp^ dxA4^^^x)sm{kiFx)sm{k2Fx). (22) 



In the last two lines, tsp{ttp) are the amplitudes of singlet (triplet) Cooper processes. 

The highly non-linear (cosine) terms in Eqs. ([l9| , ^,^l| ) signal potential instabilities of the ground state due to 
inter-subband interactions. For the dBS process [Eq.(|19[)], this instability is of the charge-density-wave (CDW) type, 
quantity (pip — ip2p being the phase of the CDW (particle-hole) condensate. If subbands are equivalent (kiF = fc2F), 
the energy is minimized by adjusting the CDW-condensate phase is such a way that the cosine takes its minimum value 
(— 1, for repulsive interactions). For non-equivalent subbands (/cif 7^ fc2F), the global minimization of the energy is 
impossible due to the position-dependent phase shift, and thus the CDW instability is suppressed. Nevertheless, if the 
energy gain due to opening of the CDW gap is large enough, the system may choose to adjust the subband densities. 
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which makes the CDW instabihty possible. Density adjustment is most hkely to occur if the cross-section of the wire 
is approximately symmetric. For example, if it is a perfect square, the second subband of transverse quantization is 
doubly degenerate. Deviations from the ideal shape lift the degeneracy but the energy splitting between the states 
remains small for small deviations. Such states are almost equally occupied, and a small difference in densities is likely 
tOphe eliminated by opening of the CDW-gap. It seems that cleaved edge quantum wires investigated by Yacoby et 
al.y satisfy this requirement. ■— ■ 

In the context of two-leg Hubbard ladders, the CDW-process of this type is known as "deconfinement"ll3: degeneracy 
of bonding and antibonding subbands implies that the amplitude of interchain tunneling, t± , is renormalized to zero 
by interactions, and electrons thus remain "confined" to their respective chains. 



The xBS process (20) contains both harmonic terms, arising from backscattering of electrons with parallel spins, 
and cosine terms, arising from backscattering with antiparallcl spins. The lattep-pontain only spin fields and thus can 
lead to the instability only in the spin channel. In the terminology of Ref.t§EJ, this instability corresponds to-jthe 
"orbital antiferromagnet phase" (OAF). The OAF instability occurs only if backscattering is sufficiently strongE^EJ. 
For a quantum wire, in which all amplitudes are given just by the corresponding Fourier components of the same 
interaction potential, this conditions means that t/(2fci?) > 2U{0) (for identical subbands), which is never the case 
for any physical f/(r). In what follows, we will not therefore consider the OAF phase. (Note that for a Hubbard 
ladder, amplitudes of various scattering processes may be determined by entirely different physics, e.g., some of them 
may result from direct electron-electron interaction and some from exchange of virtual phonons. Hence, the ratios of 
amplitudes may be arbitrary, and the OAF instability is possible, at least a priori.) 

Finally, the CS process [Eq.(|2l|)] may lead to a superconducting instability (of both singlet and triplet types), 
accompanied by opening of spin gaps in each of the subbands, in a analogy with a superconducting transition in 
higher dimensions. Quantity 9ip — 02p plays the role of the superconducting condensate phase. Inter-subband forward 
scattering [Eq. (|l^)] has an important role in developing a superconducting instability-it reduces electron repulsion 
in the relative charge-density fluctuation channel, makins-it possible for Cooper scattering to become relevant. The 
superconducting phase is also knowij_as_a "ClSO-phase"ta (meaning: one gapless charge mode and no gapless spin 
modes) or a "d-wave superconductor"Eil^t£l(indicating that the order parameter is odd upon interchanging the electrons, 
forming a Cooper pair, between subband) . This particular instability received much attention recently as one of the 
models of HTC superconductivitjtJ. We also note in passing that the idea of superconductivity in a two-band system 



has a long history, starting from the 1968 paper by FrolichEj (for a review, see Ref.H). The idea, employed in earlier 
work, is that if the masses of electrons in two subbands are significantly different, there always-even in 3D-exists 
a gapless plasmon excitation (a direct analog of Langmuir-Tonks ion sound waves in plasma), which serves as a 
mediator of effective attraction. The superconducting phase in a ID two-band system is already composed of gapless 
excitations and is not limited by the condition of different masses (although, as we will see shortly, there is no lack of 
other constraints). 

Note also that Uf (|l9|) and (21) mix charge and spin modes, and thus spoil the spin-charge separation present 



in the Hamiltonian of a single subband. 

Processes (|l^ - 21) have been written down in the literature in many different ways, so it is worth to make a 
connection to previous .prprk here. Identification of our notations with the "y-ological" ones (used by Schulz in his 
two important papersliJIlj) is as follows: tsp = gi2 + .923, Up — 923 —912, bi = (713, and his gn-process (backscattering 
with opposite spins) should be equated with our U{2knF) in Eq. ([TtI) . There is no correspondence to our amplitude 62, 
which describes exchange inter|Siibband backscattering of the type RligRl^^, Rns' Rms + {R ^ L), see Eq.(p^. That 
such a process is absent in Ref.EJ is clear from Eq.(2) of this reference. 

One more connection is made by noting that (somewhat lengthy) Eq.(pG) can be represented compactly as 



t/f = - y dx{ibi + b2){piP2 + S1S2) + [bl ~ b2)Ucl3c2 + JslJs2)}, (23) 

where p„ (Sn) is charge (spiiiL density, and jcn (jsn) is charge (spin) current in the n-th subband, using notations of 
Emery, Kivelson, and ZachaiO. 

III. SPINLESS ELECTRONS 

1. Model 

In this section we consider a "toy" model of a two-subband system of spinless electrons, which contains all interesting 
effects we want to discuss and, at the same time, allows for a rather complete analytic treatment. In this model, 
various parts of the Hamiltonian reduce to 
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Hn Hi 



0^ J dx'^ 



(24) 



(25) 



fjB _ fbs 
d ^ 2 2 



c7^ = 



(ixcos[v47r(i^i - (^92) + 2{kiF - k2F)x\\ 



(&i + b2)dxVi9xV2 + (bi - 62) fx 6*1 6*2 



(26) 
(27) 



so that 



fc 



dxcos v47r(0i — O2)] 



H = Hq + U" + + 



(28) 



(29) 



Amplitude fc plays now the role of tsp for spinless electrons. Analysis of potentially "dangerous" (in a sense of inducing 
instabilities) intersubband processes reduces to estimating the scaling dimensions of corresponding cosine operators in 
terms of the parameters of the harmonic part. In their turn, these parameters are related to the Fourier components 
of the electron-electron interaction potential. As it turns out, the latter relation is not that straightforward, and we 
will clarify this point in the next Section. 



2. Galilean invariance and Pauli principle: single- subband Luttinger liquid 



To begin with, we consider the simplest case when there is no intersubband interaction and the Hamiltonian is 
given by the sum of two single-subband Hamiltonians (p^). (As our discussion is referred now to a single subband, we 
suppress temporarily the subband index.) For a given effective ID interaction potential U {x), the Luttinger-liquid 
parameters {K and v) depend on the g = and q = 2kp Fourier components of U, as well on bare Fermi velocity vp'- 



K = ]C[U{Q)/vF, U{2kF)/vF], 
V = vfV[U{0)/vf, U{2kF)/vF] 



(30) 
(31) 



where IC{x,y) and V{x,y) are some dimensionless functions of their arguments. Relations (p0|),(|3l|) have to satisfy 
(i) the Pauli principle and (ii) Galilean invariance. The Pauli principle for spinless fermions means that for the 
case of contact interaction, i.e., when U{0) = U{2kF), the system should behave as if there is no interaction at all. 
Accordingly, K = 1 and v = vf for this case, or 



IC{x, x) 

Galilean invariance stipulates that Kv = vp, or 



V{x,x) 



K^{x,y)V{x,y) = l,yx,y. 



(32) 



(33) 



Physically, condition (|3^) comes about either by requiring thjat the shift of the ground state energy due to the motion 
of a system as a whole does not depend on the interactiorCij-ar by requesting that the dc conductivity of a uniform 
system^ does not depend on interactions (Peierls theorem)EJca. (Also, one can use the interaction-invariance of the 
persistent current in a ring threaded by the Aharonov-Bohm flux). 



° Here we consider a uniform Luttinger liquid. The role of reservoirs, to which the wire is attached to, will be discussed in 
Sec.0. 
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Conventional bosonization of the g-ology Hamiltonian (see, e.g., revie J) leads to 



K 



2ttvf + g4 - ff2 
2TrvF + 94 + 92 ' 



If-. , .g4-.g2w-, , 54 + 52^ 
w = ufW (1 + )(1 + )■ 



2'kvf 



(34) 



In terms of the Fourier components of the interaction potential, the ^-parameters are expressed as — U (0) (right- 
right and left-left amplitude) and 92 — U{Q) — U{2kp) (right-left amplitude), and Eq.(M) gives 



K 



1 



Lf(2fcf ) 



2U(0)-U(2kF) 



1 + 



V = vf 



27ZVF 

U{2kF) 

2'KVF 



1/2 



1/2 r 



1 + 



2U{0) ~ U{2kf 

2'KVF 



1/2 



(wrong) 



One can see that expressions above do not satisfy conditions (|3^), (|33|). Indeed, it follows from Eq. (wrong) that 
V ^ Vf for J7(0) = U{2kF) and that Kv ^ vf as long as U{2kF) 0. Usually, spinless Luttinger liquid model does 
not include backscattering explicitly. The rationale for such a sirnalification is that for spinless particles in ID this 
process is indistinguishable from forward scattering, see , e.g., Ref.E3. We do not find this approach satisfactory, as it 
is clear that the behavior of the system should be determined both by forward and backward amplitudes. Also, correct 
expressions for K and v should include both J7(0) and J7(2fci?), otherwise the Pauli principle cannpt be satisfied. This 
argument can also be re-phrased in terms of direct and exchange contributions to the self-energjc3. 

What did we do wrong to arrive at the Luttinger-liquid model which does not satisfy two basic physical principles? 
As one can shoaz. by using the Ward indcntities (conservation laws) for the system of interacting electrons with 
linear spectrumtZl, the problem occurs already at the level of fermions and is thus not inflicted by some subtleties 
of bosonization. Rather, it is a manifestation of an anomaly^ i.e., a violation of the conservation law caused by 
regularization, which one is forced to used in a model w ith line ar and unbound spectrumE3. 

One way to deal with this problem is to replace Eqs.( wrong) by expressions which do not follow directly from the 
original feriniO|n Hamiltonian with linear spectrum, but do satisfy all necessary criteria. This is an accord with the 
point of viewc3 that one should consider K and v as phenomenological parameters, which are renormalized from their 
bare values by irrelevant or marginal operators neglected in the course of linearization. It turns out that one can find 
exact expressions for K and w, satisfying a minimal set of requirements. 

First, we notice that the Pauli principle requires IC{x,y) to be a function of either x — y or x/y. The latter choice 
contradicts to the requirement that K, must have Taylor expansions both around x — Q and y = Q. Therefore, 



IC{x,y) ^ k{x - y), (35) 
k(0) = 1. (36) 

Then we notice that the model with forward scattering only, i.e., the original Luttinger modeS, does respect Galilean 
invariance. Therefore one can take the Luttinger model expression for K as the correct one, which means that 



JC{x,0) = 1/VT+x. 



Combining Eq. {pm with Eq. (BTj), we see that 



l^ix,y) = 1/^/1 + X - y, 
V{x,y) = ^l + x~y, 



or 



K ^ 



v — Vf 



U{0) - U{2kF) 



-1-1/2 



'kvf 



U{0) - Ui2kF) 



-,1/2 



'KVF 



(37) 

(38) 
(39) 

(40) 
(41) 



The physical meaning of Eqs. ( 40 ) , (pl|) is obvious: the effective interaction is equal to backscattering minus forward 
scattering. One can check that a Luttinger-hquid model with parameters given by Eqs. (E^, (Ell) reproduces correctly 
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results for a ID electron system, obtained without linearization but in the Hmit of weak interactions. For instance, 
the (inverse) compressibiUty of a Luttinger Hquid, parametrized by K and v from (^ , ^l| ), is given by 

1 7TV 

- = — = TTVF + U{0)-U{2kF). (42) 
X K 

As one can check, Eq. ( p2| ) coincides with the inverse compressibiHty of electrons with a quadratic spectrum obtained 
in the Hartree-Fock approximation. A perturbative (linear in U) form of Eqs.(|4^,^) has receptly been derived incil. 
It can also be read off from the tunneling exponent of a ID system with a quadratic spectrumc3. 



3. Galilean invariance and Pauli principle: coupled suhhands 

Now we allow for harmonic coupling between subbands, i.e., take into account intersubband forward [Eq.(p5|)] 
and exchange backscattering [Eq.(p7|)]. For the contact interaction, the amplitudes of these two processes coincide: 
hi — 1)2 — fo and, as a result, inter-subband interaction drops out. The Pauli principle is thus satisfied. Intersubband 
exchange backscattering does violate Galilean invariance, and the correction procedure, similar to that for a single 
subband, is necessary. We will not do it here however, because for a long-range interaction {a/d ^ 1), the violation is 
"weak": the deviation from the Galilean- invariant result is proportional to the exchange amplitudes, which are small 
compared to the direct ones. 



A. Nearly equivalent subbands 



First, we discuss the CDW-instability, which may occur if the density equilibration between subbands is energetically 
favorable. To simplify the discussion, we consider the case of a long range interaction [a/d 1), when amplitudes of 
exchange processes are small. In the leading order, fc = &i,2 = and the only "dangerous" process to be considered is 
direct intersubband backscattering. Furthermore, we assume that subbands are nearly equivalent and put vip — V2f 
and Ki = K2 but keep 5kp = kip — k2F m Eq. ( p6| ) finite. It is convenient to introduce symmetric and antisymmetric 
combinations of boson fields 



ipi ± (p2 



V2 ' V2 ' 

which correspond to fluctuations of total (-I-) and relative (-) subband charge and current. In terms of these fields 



(43) 



H 

H+ 

iJ_ = - 



dx 



v+K+id.. 



cos[v87r(y3_ + 25kFx\ 



where 



2C/(0) - U{2kF) 



1-1/2 



ttvf 



U{2kF) 



1-1/2 



ttvf 



(44) 
(45) 

(46) 



and = vf/K±. Note that K_ > 1 for U{2kF) > 0, which signals effective attraction in the (— ) channel. 



1. Collective adjustment of densities as a commensurate-incommensurate transition 



To understand how the CDW-instability works, we consider first a model situation of < 1, so that operator 
cos[-\/87r(/)_] is relevant in the RG sense. Finite Skp stops the RG-flow at scale £ ^ l/ln|5fci?|a, thus precluding the 
system from reaching its strong-coupling limit. However, this consideration does not take into account the possibility 
of a collective density readjustment between subbands. Such a readjustment may occur, if the kinetic energy loss 
fl ~ vpSkp is compensated by the gain in the potential energy due to opening of the gap in the (-) channel. In other 
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words, when the difFerence in electron densities is sufficiently small, the total energy is minimized by equating the 
densities and opening the charge gap. Obviously, such a process cannot be considered at the level of single-particle 
description of transverse quantization. Instead, one should now treat eigenstates and eigenenergies of the wire as 
being determined by a self-consistent procedure, involving both single-particle and many-body effects. 

The mechanism described above can be considered as a commensurate-incommensurate transition. The incompign.- 
surability, defined as I = ^ (J dxdx ) , where L is the length of the wire, is known to have a threshold behaviorocj: 
X ~ -y/ri^ — Q1Q{Q — flc), where flc = ^ll-nK I^np w and the expression for the gap follows from mapping on the 
exactly solvable model of a Heisenberg spin chairJlj, Acd^^ ~ (/fjs)^^^^^"^"''- As follows from the definition of the 
incommensurability, X = implies bkp = 0, i.e., equal subband densities. Therefore, the re-adjustment takes place if 
17 < r^c. Backscattering is then enabled and relevant (for K_ < 1), even for a non-zero initial value of 5kp. 

What is the physical meaning of this instability? A simple picture can be obtained in the limit of strong (both 
inter- and intra-subband) interactions, when the potential energy dominates over quantum fluctuations. In this case, 
electrons of each of the subbands form a regular lattice (Wigner crystal) . Boson fields (p„ also have periodic structures 
with period equal to y/lr (recall that a shift of -Jtt corresponds to adding one electron to the system). For fbs > 0, 
the energy of intrasubband repulsion 

fbs cos[\/87r^_] = -fbs cos[\/47r((/3i - <^2 + \/7r/2)] (47) 

is minimized by a relative phase shift of -v/tt/^ between the subbands, which corresponds to a shift of electron lattices 
by half-a-period. This is an inter-subband charge-density wave (CDW). 

V5F 



FIG. 5. An illustration of the charge density wave in two coupled subbands. A staggered configuration lowers the energy 
due to short-range repulsion, if the densities are commensurate. 



2. Competition between CDW and Cooper channels 

Let us now suppose that the density re-adjustment did occur, i.e., 5kp — 0, but Cooper scattering is also present, 
so that the Hamiltonian of the (-)-channel is 

H^ = \ [ dx{^{dx^^f + v-K_{dx0^)^ + ^cosVS^^. + -^cosVs^eX . (48) 
2 J [ K- TT^a^ ZTT^a^ J 

Which of the two instabilities-CDW or superconductivity-wins? The situation of this typjO. when cosines of both 
mutually conjugated fields (iy9„ and 9-) are present, was analyzed by Schulz and Giamarchio They found that the 
result is very sensitive not only to the value of K^, which determines the scaling dimensions of the fields, but also 
to the ratio of amplitudes, fc/fbs- As K- > 1 for repulsive U{r), it may seem that superconductivity is favored 
over CDW. The situation is not that straightforward, however. For example, consider the situation of weak and 
long-range interactions, i.e., assume that U{0),U{2kp) <C vp and a d. Because the interaction is weak, both 
processes are almost marginal, CDW being on the slightly irrelevant and superconductivity on the slightly relevant 
side. For long-range interactions, fbs ~ U{2kp) and fc ~ {a/d)'^U{0). Modehng U{x) by 



C/(x) 



_ J e^/ex, for a; < d; 
0, for X > d. 



we get U{2kF)/U{0) ~ ln{kpa)/ ln{d/a). Thus 

fc Ind/a 



fbs ^dJ \nkpa 

The RG-equation for 



< 1. (49) 



lK_^f^^-fl (50) 
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shows that K- decreases, if \fc\ < \fbs\- Even if initially K^{0) > 1, the situation with K-{1) < 1, when CDW is 
relevant, will be reached in the process of renormalization. For weak and long-ranged interactions, CDW thus wins 
over superconductivity. 

If interactions are not sufficiently weak and/or long-ranged, only a full RG solution can determine the leading 
instability. We will not analyze the general case here. 



B. Non-equivalent subbands: Renormalization Group 



Now we consider a generic situation of non-equivalent subbands, when dkp is not small enough for the density re- 
adjustment to occur. We find that a strong imbalance between Fermi- velocities of occupied subbands actually helps 
superconducting instability to develop (see case B below), despite the fact that a naive scaling dimension estimate 
does not show this. This effect follows from the next-to-leading order perturbative RG calculations, which we present 
here. 

Because Skp ^ 0, we neglect the dBS process [Eq. from the outset but keep the Cooper one [Eq. (^]. 

For long-range interaction, one can also neglect the xBS process, Eq.(p7|), whose amplitude is small for this case: 
bi oc {a/d)^. Its inclusion is straightforward (J7^ is quadratic) but does not lead to any qualitatively new results, 
while complicating the analysis significantly. The Hamiltonian then reads 



H = Hn 



(51) 



Because U'~^ contains 0-fields, it is convenient to switch from the Hamiltonian to the Lagrangian approach and to 
integrate out the (p-ficlds. The quadratic part of the resulting action is diagonalized by the following transformation 



Ml 

^i2 



cos (3 sin (3 
— sin (3 cos fi 



(52) 



where 



tan 2/3 = 



,2 „,2 ' 



Uo 



2/o \/viKiV2K2/tt 



1/2 



(53) 



and 



Ml 
^J■2 



In terms of new fields, the action is given by 



cos /3 sin f3 



sin (3 ^ cos (3 
\Jv\Kx \Jv2K2 



(54) 



S 



dxdr 



E 



Rn{ — idr 



■ COS 



471(6*1 



(55) 



where 



4? 



(56) 



are the velocities of new collective modes and i?„ = l/(u„/i^) are the new stiffness coefficients. 

We are now ready to perform the momentum-shell RG, i.e., to expand perturbatively in coupling constant fc and 
integrate out high-energy modes with 2-momentum k within a thin strip A — dA < k < A (dA/A <C 1). The first-order 
contributions renormalize /c, whereas the second-order ones renormalize stiffnesses i?,|_.The main difference between 
our case and the conventional RG-treatment of the sine-Gordon action (see, e.g., Ref.cZl) is that the /^-contribution 
produces (among others) mixed gradient terms of the type d^9id^92 {v = t,x), which are absent in bare action ([s^). 
To eliminate these terms, we transform fields one more time : 



= 1 + ^ 



cosh/3' t'sinh/3' 
fsinh/3' cosh/3' 



(57) 
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where 

U2R2 , dA / fcY 1 /t-QN 
t = — ^ and /3 = — — — — (58) 

UlKl A \ TT J U1U2K1K2 

are chosen in such a way that the coefficients in front of mixed gradient terms vanish. When written in terms of 
the action is brought into its original form but with renormahzed parameters. Resuhing RG equations read 

-— = - P ^ ( + 2 

dlRi ^ Ri\Rl{l+-i^) Rl{l + -i^) R1R2)' ^ ' 

d I _ pi ( 4 27^ 1 X 

dZi?2 i?2Vi??(l + 7') i?i(l + 7') ^i^27^' 



where / — (l/7r)/cy ^ + U2 ^ is the dimensionless couphng constant and 7 = ui/u2- To the /^-accuracy, ah terms 

multiplying on the right-hand-side of the first three equations above have to be treated as constants determined 
by the initial conditions. The system of RG-equations has an obvious integral of motion 

c = - + ^ - /^ (63) 

Cl C2 

where 

x=l-l/Ri, y=l-l/R2, (64) 

and C1.2 are the coefficients in front of in Eqs. (|59| , |60| ), respectively. Note also that x — {01/02) y + p (we denote 
X = x{l), y — y{l), whereas initial values are denoted by sub-index 0, i.e., x(0) = xq, etc.). Constants of motion C 
and p are determined by initial conditions. 

The flow described by ( p9[]6^ ) is quite similar to that of a canonical Kosterlitz-Thouless system: x and y increase 
with / regardless of its sign. If xq, yo > 0, f grows unrestrictedly, flowing into the strong-coupling regime with a 
gap in the 9i — O2 channel. Such initial conditions correspond to i?i,2 > 1, i.e., to the attractive interaction in the 
0„-channels. It is worth emphasizing here that due to the presence of inter-subband forward scattering, such effective 
attraction may arise in a purely repulsive system, as we shall demonstrate shortly {case A below). Another relevant 
limit is represented by the "repulsive" case (case B), where initially ccq < 0, j/o < 0. For a strong repulsion (i?i,2 ^ 1), 
/ quickly renormalizes to zero and the resulting phase is a two-subband Luttinger liquid. However, there is a region 
of anomalously small XQ^ya ^ /o (which requires strong inter-subband scattering), where Cooper scattering may still 
be important. One finds that if 

(ci + C2)/^ > {xo + yof (65) 

The Cooper process wins over repulsion and initially negative variables x, y change sign during renormalization. In 
this case /c decreases with /, initially but then passes through a minimum and flows finally into strong-coupling regime 
/ > 1. Equation ( |65| ) is a condition for the development of superconducting fluctuations in the system with purely 
repulsive interactions. 

Now we again apply our analysis to a wire with weak and long-range interactions. Two limits are possible. 
Case A: H^cdw /kp viF — V2F < Uo- 

The first inequality allows one to neglect direct backscattering, which leads to inter-subband CDW, whereas the 
second one allows to consider subbands as "nearly equivalent". Denoting vp = fiF, Svp = vip — V2f, and the 
2fci;'-component of the interaction potential in the n-th subband by C/2feF ' finds 



(1) , rr(2) r.. ttW 

< 2, (66) 



1 ,1 U^;/^+UZ SvpU,, 



i?l i?2 2nVp Vp ^TTVp 
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which corresponds to effective attraction. Thus, thanks to inter-subband forward scattering, the Cooper process is 
relevant in the system with purely repulsive interaction. Note that a small velocity imbalance Svp > enhances the 
relevance of Cooper scattering . (As our second subband is chosen to have a higher energy, Svp is always positive.) 

Case B: vip — V2f S> U{0). Van Hove singularity. 

In this limit an analytic solution is also possible. Generally, one finds that l/i?i + l/i?2 > 2, which corresponds 
to effective repulsion. Neither CDW nor superconducting instability can develop, and the resulting phase is a two- 
subband Luttinger liquid. This is not true, however, in the limit of a strong velocity imbalance, when V2f/vif 1, 
i.e., when the second subband just opens for conduction. Then u„ w u„ and 7 = ui/u2 ^ /v2f S> 1. Hence ci ~ 7 
and it follows from ( |65|) that Cooper process wins over repulsion, if fc > [/(0)/y^. For long-range interactions, the 
last inequality reduces to 

ViF ^d/ 

The physics of this scenario is well-known - interactions are enhanced due to the large value of the density of states 
(cx 1/v2f) in tbe upper subband (Van Hove singularity). We should warn here that our calculations do not describe the 
very onset of conduction in the upper subband, because its proper description requires accounting for the nonlinearity 
of the electron spectrum, which is beyond our bosonization analysis. However, such a calculation was performed by 
Balents and Fisherlla, who analyzed the case of a contact interaction. They found that superconducting fluctuations 
are indeed enhanced in this limit. 

We thus see that superconducting fluctuations do have a good chance to overcome the electron-electron repulsion 
and drive the system into a strong-coupling phase with the gap in the spectrum of relative current fluctuations, 9i — 02- 

There is another important feature of the RG-flow described by Eqs. ( p9[]6^ ): the interaction tries to equilibrate 
densities in the subbands. This is seen from equation (|6l|): dj/dl is proportional to 1 — 7^, which makes 7 = 1 a 
stable fixed point. 7 tends to increase, if initially 70 < 1, and it tends to decrease, if 70 > 1. 



IV. ELECTRONS WITH SPINS 



Guided by the results of the previous Section, we now comment briefly on what happens if spin is included. As 
should be clear from the complexity of Eqs.(|l6[]2l|), this question has no simple answer. For a quantum wire with 
< U{2kF) < U{0), possible phases are again (i) Luttinger liquid, (ii) inter-subband CDW, and (iii) Cooper phase 
(superconductor). As with spinless electrons, subbands must be nearly equivalent in order for the CDW phase to 
occur, whereas the Cooper phase needs effective attraction in the relative charge-density excitation channel. When 
neither of these conditions is met, a two-subband Luttinger liquid is realized. At different degrees of generality, 
renorraaJizaJdoD group analysis of the model defined by Eqs.(p^ ^ ) has been performed in the past and we refer to 
papersOEjiljO for a detailed description. 



1. Long-range interactions 



If the interaction is long-range and weak, a considerable simplification occurs. In this case, amplitudes of forward 
intra- and inter-subband processes are the same [see discussion after Eq.(p^] and a simple perturbative estimate of 
scaling dimension Sc of the Cooper process (^ ) is possible. 

For weak interactions, Knp = 1 — {2U{0) — U^^) l2'n:VnF, and /o = U{Q). The S'C/(2)-invariance requires that 
Kna — 1- One finds 

5^^2~^^Hi±^-^^<2. (68) 
Attvf vf 4:Ttvf 

Observe that this result coincides with Eq.(^6|) upon replacing U{0) 2U{Q). Thus, Cooper scattering is relevant for 
repulsive long-range interactions (and assuming also that < Svp ^ Uq). 

However, if Kp^ is sufficiently close to its non-interacting value, i.e., to unity, backscattering is strong (/^s 3> ts^^^t^ 



and 5kF is small, the CDW-channel can take over the Cooper one, similarly to the scenario described in Sec. Ill A 2 
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2. Electron ladder 



There is an important question where an RG consideration is very helpful. Suppose that condition ( pq ) is satisfied 
and thus tunneling of fluctuational Cooper pairs is relevant. What happens to spin excitations? To answer this 
question, we relax the SU (2)-invariance condition and perform the RG calculation for two nearly equivalent subbands 
so that Vnu = Vu, Knv = K,y, where v = p, a. Nevertheless, we assume that Skp is still finite and neglect direct 



backscattering (19), similarly to Subsection [II A . The problem then becomes identical to that of two coupled equiva- 
lent chains coupled by the interaction ("electron ladder"). Also, for the sake of simplicity, we consider only the singlet 
channel of Cooper scattering. Due to enhanced symmetry, the total current fiuctuation mode 9p+ — {9ip + 92p)/V^ 
decouples from the rest and is described by a harmonic action with 

Kp+ = [K-^ + 2/o/^«f] . (69) 

Relative current fluctuations 6'p_ ~ {9ip — 9ip)l\f2 are described by the sine-Gordon theory [t^p term in (|2l])] with 

K-'^ = _ 2/o/(^z;f)] and Vp- = Vp{\ - 2foKp/TTVFY''' ■ (70) 

In addition, we have to keep track of the spin-density sector, which contain the cosine term corresponding to intra- 
subband backscattering (|l^). For convenience, we denote the amplitude of this term by its initial value being 
5(t(0) = U{2kF)- After tedious but straightforward calculations we arrive at the following system of RG equations: 

^^g = 2{l-K,)g-P, (71) 
|t-=(2-X.-^-^)t-, (72) 

^(l-^.)-^(.9' + ?), (73) 

where g = ga/i'KVF) and t = tsp/ {t^vf)- Let us recall what happens in the absence of Cooper tunneling first and 
set tsp — everywhere in this system. In the weak-coupling limit, Ka = (1 — 5)^^^^ « 1 -f g/2 and ( [7l| ) becomes 
dg/M = —g^, which gives gt = go/(l + 3o^)- For repulsive interactions (go > 0), 5 oc ^ as ^ — > cx): intra- 
subband backscattering is marginally irrelevant. Observe now that when the Cooper process is present and relevant, 
i.e., when t increases, the flow of g is modified: the P-term on the right-hand-side of Eq. ( [7l| ) changes the sign of 
dg/d£. Hence gi is bound to become negative in the process of renormalization and grows unrestrictedly in its absolute 
value. Intra-subband spin backscattering is thus driven relevant by singlet-pair tunneling, which results in pinning of 
ipa in Eq.(p7|) and opening of the spin gap. Thus, similar to the true superconducting state in higher dimensions, the 
Cooper phase is characterized by gaps in both charge- and spin-channels. The only massless excitations are those of 
the total charge channel. This phenomenon is not restricted to the degepecate electron ladder but rather is a generic 
feature of the system of coupled subbands and/or chains, see, e.g., RefsJl3'tZl. 



V. CONDUCTANCE 



Having realized the importance of inter-subband interactions, we now proceed with the analysis of its effect on 
observable properties of quantum wires. The first property we consider is conductance G. 



A. No disorder 



Our results for the conductance of a clean wire can be understood from the following simple considejLations. The 
dc conductance of a single-subband wire is equal to 2e^ /h regardless of the interactions in the wireC3E£l. Consider 
now a wire with several subbands occupied. Those interband interactions, which do not open gaps, lead only to a 
renormalization of Luttinger- liquid parameters. As these parameters do not enter the final result for G, the conduc- 
tance remains quantized in units of 2e^//i per occupied subband. Other processes, such as direct backscattering and 
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Cooper scattering, may open gaps in channels of relative charge fluctuations as well as in spin channels. Neither of 
these gaps, however, affects the center- of-mass of the electron fluid through the wire, which continues to move along 
the wire unrestrictedly. Therefore, G remains unrenormalized by this type of interactions as well. 

Now we demonstrate the proof of the statements made above. Consider the case of a superconducting instability, 
when the cosine term of Cooper scattering in Eq. ([l^) is relevant and the 6'p_-field is thus gapped. Gaussian fluctuations 
of the gapped field can be described by expanding the relevant cosine term around its minimum value: 

(Afc/n^a^) cos[V2^(0ip - 02p)] » const + m^Oip - 92pf , (75) 

where m is the mass of the field. The strong-coupling (superconducting) phase corresponds to m 7^ 0, whereas in the 
Luttinger-liquid phase m — 0. Expanding cos[^/2nipna] around their minima as well, we find that at the Gaussian 
level charge and spin modes decouple again. Note though that now these are massive modes. 

As spin excitations do not affect charge transport, we concentrate on the charge sector of the theory, whose 
Hamiltonian is given by the sum of Eqs. ([l6| , |l8| ) and (|75|). To simplify notations, we suppress index p in this section, 
so that 9^, etc. Using VnKn = VnP, we write the charge Hamiltonian as 

+^ [{d^V+f - {d,W-f] + {viF'V2F)dJ+dJ- + _ d^^+d^^^^+m'eiy (76) 
The total charge current is given by j — e.yj2/iT^^ dtifn — e{2/ ^/Tr)dt(p+. The conductivity 

a{x, to)^- f dt'{[j{x, t')j(0, 0)])e'"'' = {2e/V^f{tLu)G++{x, lu) (77) 

is related to the retarded Green's function G+-^-{x,t) = —il^{t){[(p-^-{x,t),lp^ (0,0)]), which is coupled to another 
Green's function G \-{x,t) = —iQ(t){[ip^{x,t),(p+{Q,0)]) by the following equations of motion 



+ 



\2Vi^i K2J TT J + 2\Ki K2) ^ ' 

In the massless limit (m = 0), this system of equations is solved readily. In order to model the effect of non-interacting 
electron reservoirs, which the wire is attached to, we assume-ithat Kn, Vn vary with x adiabatically and approach their 
non-interacting values of Kn = 1, w„ = vp for x — > ±cx)E§"l2I. In the zero-frequency limit, the solution is particularly 

simple: G'++(x, w — > 0) = l/2iw, G y{x, w = 0) = 0. As a result, conductivity is x-independent, and the conductance 

is simply G = 2 x 2e^ /h. 

In order to see the effect of the gap, we consider first the case of equivalent subbands ("electron ladder"), introduced 



in Sec. IV 2. One observes immediately that conditions vi — V2, Ki = K2 lead to complete decoupling of equations 

for G++ and G |-. As a result, the total charge mode (p+ is not affected by the gap. Taking the boundary condition 

for K and into account gives again the universal result G = Ae^/h. The result for the "electron ladder" thus gives 
us a hint that G remains at its universal value despite the presence of the gap in the relative charge channel. In order 
to prove this statement in the general case, we neglect for the moment the boundary conditions for Kn and Vn-, i.e., 
consider a uniform wire with two coupled subbands. System (|7^ ) is then solved by Fourier transformation. The key 
feature of the result for G++((;, w) is that it still has a pole corresponding to a massless mode w oc q, despite the 
presence of the massive term. The conductivity becomes 

^^(^,g) 



'^il^'^)- r.2 ,-,2„2V,„2 j_ ,,,2 _,,2„2V (^9) 
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where w, u are some renormalized velocities, ujm is some energy proportional to m? , and q) is a smooth function 
of its arguments. As a result, a{q,u}) = G6{q) in the limit w ^ 0, where G has a meaning of the conductance. 
Because we neglected the boundary conditions corresponding to presence of non-interacting leads in this calculation, 
G depends on all interaction parameters - Vn,Kn, and m^-and is of course different from Ae^/h. It is clear though 
that once the boundary conditions are restored, this non-universal value is replaced by the universal factor of Ae^ /h. 
The only other possibility is G = 0, which, however, is ruled out by the fact that G++{q,uj) has a masslcss pole. 

We thus conclude that the conductance of a clean wire remains at the universal quantized value irrespective of 
whether the relative charge mode is gapped or not. The case of a CDW instability can be treated in a similar manner. 



B. Disordered wire 



A disordered two-band system in th(Hpresence of interaction-induced instabilities was considered by Orignac and 
Giamarch£3 and by Egger and Gogolinc2l. Our discussion of a disordered two-subband wire follows largely these two 
papers. 

Results of the subsequent analysis can be summarized as follows. If Cooper scattering opens a gap, the system 
does not become a real superconductor: a single weak impurity splits eventually the wire into two disconnected halves 
at low enough energies, and even a weak random potential leads to localization of electrons, similar to the case of a 
gapless Luttinger liquid. Nevertheless, effects of disorder are less pronounced than for a gapless Luttinger liquid. On 
the contrary, the CDW-state is more sensitive to disorder than a gapless Luttinger liquid. 



1. Spinless electrons 

We begin by considering a single impurity described as a potential perturbation w{x,r±). The impurity causes 
backscattering of electrons within the occupied subbands, as well as inter-subband backscattering. The amplitudes of 
corresponding processes are given by 



Wn{2knF) — J dxdrj^w{x,r±)(l)^{rj^)cos{2knFx),n—l,2; (80) 
Winter = dxdrj^w{x,rj_)(t)i{rj_)(j)2{ri_)cos[{kiF + k2F)x] . (81) 



If w varies slowly across the wire, then Winter ^ W„ due to the orthogonality of transverse wavefunctions, and we 
consider intra-subband backscattering first. The bosonizcd form of intra-subband backscattering is 

Wn{2kr,F) , rr 1 W„(2A:„f) 



= cos[V4^^„] = ^lllliif^ cos[V2^(^+ ± ^_)], n = 1, 2; (82) 

so the total backscattering operator is given by 

W2{2k2F) - Wi{2kiF) . 



Wintra = > , Wintra = ^ ' ^ Sm[V 2tT ip +] Sm[V2TT(p^] 

„tt2 

, W2(2fc2F)+Wi(2fciF) 



na 



cos[v27r(y9+] cos[v27r(^_]. (83) 



Note that Wintra is local in space and thus cannot change the RG-flows of bulk parameters of the wire. Depending on 
these bulk parameters, however, Wintra will either grow, splitting eventually the wire into two disconnected pieces, or 
decay, in which case the impurity effectively disappears, 
(i) Cooper phase. 

In the Cooper phase, 6- is gapped, hence (p- fluctuates strongly. [This follows from the fact that 6- and ip- are 
canonically conjugated fields, see Sec. On the first sight, it may seem that these strong fluctuations render Wintra 
to zero. To see that it is not so, consider a second-order impurity contribution, e.g.. 



/ Wi(2fcif) 



^ y"dTy"dT'(e'^'^-(^)e-*^'^-(^'))cos[%/2^(p+(r)]cos[V2^^+(r')] (84) 
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As it will be explained in more details in Sec. VI, the correlator of (/5_-fields in the Cooper phase decays exponentially, 
i.e., i^Q^^^'f-i'^)er'^"^'P-^^''i) ~ Q-'^sc\r-T'\ ^ where Asc is the Cooper gap in the 0_-channel. As a result, the double 
integration over t, t' in Eq. (|4|) is effectively contracted into a single one, the result of integration being 

Ag^{W{2kiF)/naf J dr cos[VS^(p+{t)]. (85) 



The mechanism of generating higher order impurity backscattering was discovered incjiij. (vStt under the epsine 
indicates that this is a two-particle backscattering process.) Following the RG-calculations of Kane and FisherEI, one 
finds that impurity backscattering, generated in this way, becomes relevant for < 1/2. Note that this requires 
rather strong electron repulsion. For weaker repulsion, backscattering amplitude scales to zero and the wire retains a 
universal conductance of 2 x e'^/h. Without superconducting correlations, i.e., when 0_ is not gapped, an impurity is 
effective for K < 1. Thus the Cooper phase weakens but does not eliminate impurity scattering, 
(ii) CD W phase. 

In this phase, Skp = and ip- is pinned by the bulk nonlinear term cos[V87riy9_], so that (p- acquires average value 
((/?_) = y^Tr/S. Allowing for fluctuations around the average value, we substitute = ((^-) + Sip- into (p3|), which 
gives 

w W2{2kF)~Wii2kF) . . /T^ 1 r /^A 1 W2{2kF) + Wi{2kF) /— , . . /.^x i ^««^ 

M^intra = ^ ^ ' sm[v27ri^+] cos[v27rdi^_] ^ ^ cos[v27ri^+] sm[v27rdi^_] (86) 

na na 

Observing that for small fluctuations one can replace cos[\/27r(5iy9_] 1, we see that the first term in ( |86| ) gives the 
strongest contribution to backscattering, which is relevant already for A% < 2. The second term requires more work. 
To the second order in the amplitude of this term, an expression similar to ( p^ ) is generated, but now it involves the 
following average 

(sin[V2^,5(p_(T)]sin[V2^<5(^_(r')]) ~ sinh[Xo(AcDw|r - t'|)], (87) 

where Acdw is the CDW gap, Aq is the ultraviolet energy cutoff, and Kq{x) is the modified Bessel function, 
[Kq{x) I \fx for x^\\. This result is due to the fact that in the massive phase 

^^^V2^SvAr)^±^V2^5^-ir')^ _ {AcDw/Ao) T Ko{AcDw\r ~ t' \)] . (88) 

As a result, correlation functions of sines and cosines of massive fields are different: the first ones decay exponentially 
with distance, whereas the second ones reach constant values. Thus the double integral over r, r' in (^) can be reduced 
to the single one again, and, similarly to the Cooper-phase case, two-particle impurity backscattering, relevant for 
K+ < 1/2, is generated. 

It is also straightforward to analyze the effect of inter-subband impurity scattering 

Winter = cos[v27ri^+] cos[v 27r6'_] . (89) 



Similarly to Eq.(84), wc have to average over the strongly fluctuating 6'_-field, which generates again the two-particle 
backscattering term ~ cos[\/87r(p+], relevant for < 1/2. 

Hence, the perturbative correction to the conductance of a CDW wire behaves as 

- SGcDW (X wh^'^-^ + {wyAcDwfe*''+-^, (90) 
where e = max{T, bias} and where we have also indicated the o rder of the imp urity potential. Please note that 



the exponent of the weak-link counterpart of (|90|), derived in Fa. (140) of Section Vll, is not related to the leading 
exponent — 2 in Eq. (|90| ) by the conventional duality relationD. We conjecture that this violation of duality signals 
phase transition separating regimes of weak and strong tunneling. 

For the Cooper-phase case, cos['\/27r0_] is replaced by sin[V27r(50_]. Hence, Winter also generates the effective two- 
particle term in the second order of perturbation theory (cf. the analysis of the second term in (p6|)). The correction 
to the conductance is given by 

-^Gscoc (z«VA5c)' (91) 
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To summarize, the Cooper phase is insensitive to a single impurity as long as > 1/2, whereas the CDW one is 
stable only for if+ > 2. 

We now turn to the case of a weak random potential produced by many impurities. To establish the boundary 
between delocalized and localized regimes for the cascjof weak disorder, it suffices to replace e ^ 1/L m Eqs. (90 9l| ) 



and to multiply SG by the total number of impuritieiil, proportional to L. Depending on whether SG increases or 
decreases with L, the wire is in the localized or delocalized phase. By doing so, one concludes that a wire is localized 
for K+ < 3, if it is in the CDW phase, and for K+ < 3/4, if it is in the Cooper phase. For comparison, a (spinless) 
single-subband Luttinger liquid is localized for K < 3/2. 



2. Electrons with spins 
The bosonized form of impurity backscattering is given by 

''^'^^ " 2^ cos[V^((^p+ ± ipp.)] 

X cos[V7r(v3<T+ ± ipa- )] ; (92) 

Winter (fclF + fc2F) , ^, , /) M 

Winter = Z C0s[V7r(<y5p+ + Op-)\ 

X cos[V^{9^- + ipa+)], (93) 

where +(— ) in the argument of cosines refers to the 1st (2nd) subband and all operators are evaluated at the position 
of the impurity. 

In the Cooper phase, the Op-- and i^cr±-rnodes are gapped, whereas the conjugated modes, i.e., (pp- and 9a-± , 
exhibit strong fluctuations. Integrating out (pp_ and ipa± , we find 



Wintra ^ A^^ [W(2fci^2 f)] C0s[V4^(^p+] , (94) 

which is relevant for Kp^ < 1. The same is true for Winter > where strong fluctuations of O^^^ produce a similar 
operator. The correction to the conductance behaves as 

-6Gsc^{w/Ascfe^''-+-^, (95) 

i.e., as if we were dealing with a single-channel Luttinger liquid, characterized by parameter Kp^, subject to an 
effectively reduced impurity potential. Weak random potential leads to localization for Kp+ < 3/2. .. . 

In the CDW-state, the situation is different. In this case, the fp-- and (pa-±-modes are gappedoEfl, whereas the 
Op-- and Oa±- modes fluctuate strongly. As a result, intersubband backscattering is renormalized into cos[-\/47r(pp-(-], 
as in the Cooper phase, but intra-subband one remains unchanged and is determined by the dynamics of the only 
gapless (fip^ mode: 

Wintra (X C0s[V7ri^p+]. (96) 

Pinning of ipp+ at the impurity site leads to the suppression of the conductance. Wintra is relevant, i.e., an impurity 
eventually splits the wire into two disconnected halves for Kp+ < 4. The correction to the conductance behaves as 

-SGcDW ocw'^e^-^. (97) 

This is to be contrasted with the case of a gapless Luttinger liquid, when the impurity is relevant only for repulsive 
interactions {K < 1). This reflects the fact that a real (gapped) charge-density- wave is pinned stronger than the 
fluctuating one (Luttinger liquid). Finally, weak random potential localizes the CDW- wire for Kp+ < 6. It is worth 
pointing out here that such large values of critical Kp+, separating localized and extended regimes, imply strcmg 
effective attraction between charge fluctuations in the p+ channel. It might well be that a charge segregationE3, 
instead of the CDW instability, will take place for such a strong attraction. 

Thus, for electrons with spin, the Cooper phase is more stable to impurities than the CDW one, similar to a spinless 
case, but both are unstable in the physically relevant region of Kp^ < 1. 
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VI. SINGLE-PARTICLE DENSITY OF STATES 



We now turn to the discussion of tunneling into a quantum wire. To the leading order in barrier transparency T, 
the differential tunneling conductance is 

G{V) = \T\^PcP{eV), (98) 

where pc is the density of states (DOS) in the contact (which we assume to be energy-independent), p{e) is the DOS 
of the wire, and V is the applied voltage. When the wire is in the gapless Luttinger-liquid phase, p{e) oc \€\^ . This 
behavior has recently been observed in tunneling into carbon nanotubesEL Tunneling into the edge of a fractional 
quantum Hall system also exhibits a power-law current-voltage dependences, which might be an indication of a chiral 
Luttinger-liquid state at the edge. Here, however, the situation is .pot that straightforward, and other explanations, 
different from a chiral Luttinger liquid, have also been suggestecOo. 

Suppose now that a two-subband quantum wire is in one of the possible gapped phases, i.e., CDW or Cooper 
phase. The goal of this Section is to analyze what would a tunneling experiment show in this case. The answer 
turns out to depend crucially on the geometry of the experiment. If the tunneling contact probes the interior of the 
wire, the gapped behavior is predicted: G{V) = for eV < A, where A is the appropriate energy gap. For eV > A 
the behavior is non-universal: the threshold behavior of p is determined by gapless charge and spin modes. More 
surprisingly, tunneling into the end of the CDW-wire exhibit a gapless behavior, similar to the Luttinger-liquid case. 
The tunneling exponent though is different from that for the gapless phase. 



A. Tunneling preliminaries 

The local single-particle (or tunneling) density of states is given by 

p{uj,x) ^ --Im{Gret{uJ,x)}, (99) 

TT 

where Gret{^^, x) is a Fourier transform of retarded Green's function Gret{t, x) = —iQ{t)'E,s{{'^s{t, x), 5'J(0, x)}). Rep- 
resenting the electron of the n-th subband as a sum of right- and left-movers and accounting for the orthogonality 

of transverse wavefunctions, the Green's fmiction becomes Gret{t,x) — G^^f{t,x) + G^^f(t,x) , where the sum- 
mation is over all occupied subbands and G'^^{t,x) — —iQ{t)'Es{{Nnsit,x),Nlg{0,x)}), N — L,R, see (^jl4). The 
contribution of the off-diagonal terms ~ {RnL\^ is less singular and is thus neglected. p(w, x) is a sum of contributions 
from right and left movers of all occupied subbands. To find p{lu^x), it is convenient to calculate first the Matsubara 
Green's function 

G«(t,x) = -{TrRs{T,x)Rl{Q,0)), (100) 

and then to make the analytic continuation to real frequencies. Left- and right-moving fermions give identical contri- 
butions to p, thus the result obtained from ( |100| ) is simply multiplied by a factor of two at the end. 

The key feature of gapped phases in multisubband ID systems is the co-existence of gapped and gapless modes, which 
also makes the calculations to be slightly less trivial. The single-particle Green's function under these circumstances 
has recently been considered in Ref.EScZl, and our analysis follows largerly these two papers. 



B. Warm-up: DOS of a half-filled Hubbard chain 



To warm up, we consider the simplest system in which gapped and gapless modes co-exist ~ a single-band Hubbard 
chain at half-filling. In the context of nanostructure^hysics, such a system is produced by imposing an artificial 
periodic potential of period qq over a quantum wireEj. At half-filling, the Fermi-momentum kp = 7r/2ao is com- 
mensurate with the reciprocal lattice spacing, which gives rise to Umklapp scattering. An Umklapp process occurs 
as simultaneous backscattering of two right- or left-moving electrons, the total momentum transferred to the lattice 
being ±4 x 7r/2ao — ±27r/ao. This process is responsible for opening of the (Mott-Hubbard) gap in the charge- 
excitation spectrum. On the other hand, spin excitations remain gapless, and are described by the S'[/(2)-invariant 
Luttinger-liquid Hamiltonian (K^ = 1)- 

The corresponding Hamiltonian of the charge sector is 
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H=\j dx{vpKp{dJpf + ^{d^Vpf + gcos[V%^ipp]}, (101) 



where the last (cosme) term represents Umklapp scattering. Substituting Eq. (|lj) into Eq. (100), one finds that 
factorizes into a product of spin and charge parts 

G«(r,x) = -f|^F,(r,x)Fp(r,a:); 

= {eMi\[\{M^) - exp[-iy|(e.(l) - ^.(0))]).; v = p,a, (102) 

where shorthand notations 1 = (r, x), = (0,0) have been used. is gapless, whereas Fp contains massive fields. 

1. Bosonic calculation 



To calculate Fp in the bosonic language, we adopt the semiclassical approximation, i.e., expand the cosine in ( |101| ) 
around its minimum to the second order in fluctuations. This is equivalent to replacing g cos[-\/87riy9p] — > TO^(^p, which 
defines the mass = ing. Now the averaging is straightforward: 



Fp(r, x) = exp ■ 



/ (0(1 - cos[fc . z\){Go{k) + G,(fc) + 2»G.(fc) ^ I , (103) 



where k= {ko,ki) = (u!n,q), z= {t,x) and 



vK 

G^{uJn, q) = {Tr(pp(pp)u,^q = — (104) 



The Green's function of 0p-fields can be written as 



V 1 u)'^ ifi^ 

Ge{uJn,q) = {Tr9p9p)^^q ^ 17\ 2 2 , 2 + ^ 2 2 i 2\^ 2 2 2 i - 2\ ) = '^'o^'^n, q) + Gf\uJn, q), (105) 



where m = m vK. The first term in Eq. (105) is just a free Green's function, whereas the second one is present only 
in the strong-coupling phase and contains a infrared divergence at g — > 0. This divergence is often explained by the 
"uncertainty principle" : in a gapped phase, a position-like field {(fp) acquires an average value, hence its canonical 
conjugate, momentum-like field {6p) fiuctuates strongly, hence the average {OpOp) diverges. Let us analyze the Fourier 
transform of Gg in more details, and define 



Ii^,r)^ I ^(l_cos[^.z1)Gf (c.,g). (106) 



Changing to polar coordinates cjn ~ k cos <j),q = k sin (f>,T — z cos "djX — z sin i9, we get 

(27r 2 7 '^sin2 ' ' Jq y + z^ cos^ {4> - d) 



The integral over cj) diverges at = 0, which, if taken literally, means that I = oo and thus Fp{T,x = 0) = for any 
finite T. However, this divergence is absent at •& = tt/2, which corresponds to t = 0,z = x. Let us therefore continue 
the calculation at this special point. Despite the cancellation of the infrared divergence, the integral is still controlled 
by the region of small (p: cj) ^ l/{fhx) <C 1. Expanding sine/) ^ (p and extending the limits of angular integration to 
±00, we find that I{z,d = tt/2) = rhx/A + 0{l/{rhx)) for fhx ^ 1. Collecting regular contributions from G,^ and 

G^g\ we find that the equal-time exponential correlator of 0p-fields is given by 



l^^iaOpixfi) g-mep(o,o)^ = exp 



a2 



ln(a; /a ) - — "^a; 



AttK ' ' ' 4K 



(108) 



An important feature here is that the expected exponential decay of this correlator is modified by the 
power-law prefactor, given by the usual Luttinger-liquid correlator. Symbohcally, (e^^^p^^^'^^ e~'"^''^'''°^)m^o = 
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I^ia9p{x,0) g- 



■ia9p(0,0)\ 



/m=o X e"™^. One should be careful in using Eq.( |l08| ): Luttinger-liquid parameter K, which 
appears here, should in fact be understood as the strong-coupling fixed-point value, K* , which is often unknown. 

Fortunately, the fixed-point value of K is known for a half-filled Hubbard chain: K* = 1/2. For fhx ^ 1 we then 
obtain 



0p(O)\ 



-!■ const, 

exp[— 7rma::/4]/\/x, 



(109) 



and the full Green's function behaves as 



G^(0,x) cx 



exp[— Trma; 



/4] 



(110) 



in agreement with Ref.Ell. The second x^^/^-factor in Eq. ( |110| ) is due to gapless spin excitations. 

So far, all calculations have been straightforward. Now we would like to argue that the infrared divergence of 
I{z, d ^ 0) is an artifact of the semiclassical approximation, which ignores degeneracy of cos[\/87r(p] with respect to 
a uniform shift (p (p + y^N with integer N. The proper theory of both massive and massless phases should be 
Lorentz-invariant. We thus propose that the correct result, valid for any z = \fx^ 



w^T^, is given by Eq.(108) where 



X is replaced by Euclidian distance z: x 



J^iaeix.T) g- 



9(0,0)\ _ 



z. Therefore, 
ln( 



exp 



^ttK 



X^ + v'^T^ 



^rhy/x^ 



(111) 



Similar arguments in favor of such replacement were given by VoitH. 

We now use Eq. (Ill) to evaluate Eq. (103), and find [compare with (|109|)] 



a 



FpiT-,0) ~ J;^exp[-7rmw|r|/4]. 



(112) 



The spin sector average is non-zero and un ivers al (thanks to = 1), therefore Ffj{T, 0) ~ ^Ja/v\T\. The correspond- 
ing DOS will be calculated later, see Eq.(123). The correctness of the procedure described above is verified in the 
next Section. 



2. Re-fermionization. 



To check that our proposition makes sense, we now switch gears and |d£rive Eq. (112) in a completely different 
way. To this end, we use the Luther- Emery refermionization procedureEj, which works for K = 1/2, i.e, at the 
fixed point of a half-filled Hubbard chain. This procedure begins with an innocuous looking transformation ipp = 
Op — \/29, which changes Umklapp scattering in Eq.( |lOl| ) into backscattering of some auxiliary particles 
(solitons): cos(\/87r(y3p) = cos(\/47r(/3). Right- and left-going solitons arc defined by 

^± = ^^e^p{±iV^iipTO)} . (113) 



In terms of new bosons, the original fermion operator (|13|) becomes 



Rs = ^ — exp[«sW -{ipa - Oa)] exp(«— (/?) e-xp{-i^e). (114) 
v27ra V ^ ^ 

It can also be written in terms of solitons 

Rs = e*^/« eMis^j\{^a ~ 9,)] exp(-z^(^) (115) 

where v? is expressed in terms of soliton density fluctuations as {\ / =: + ''P-ip- '■■ The usefulness of 

these formal manipulations is based on the fact that Hamiltonian ( |10lD is quadratic in massive solitons V'±, and the 
mass (or the gap A) in their spectrum is determined by coupling constant g: A — g/{2na). Due to the presense of 
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the gap, charge density fluctuations are suppressed, which results in the suppression of fluctuations of ip. Therefore, 
at energies below the gap the phase factor exp(— 1-^/711^9/2) in Eq. (11^) can be replaced by its average value. The 
charge part of the Green's function [Fp] coincides with the Green's function of massive fermions 



dujndq 
72^' 



vq 



(116) 



dcOn 



= = dr 

A2 V 



dtOn 



A 



sgn{T)Ki{A \t\ 



(117) 



At X = 0, 

Fp(T,0)~- 

Asymptotically, Fp(T, 0) ~ (A/u-\/A|r|)e^'^l'^l , in agreement with our earlier proposition ( |112| ). 

1/2 

It follows from Eqs.(116)-(117) that upon continuing to real frequencies i^^'^*(ijj — icj„) = —{lo/2v) (A^ — cj^) 
Hence, the DOS of massive fermions is given by 



cos(cj„r) 



(118) 



Up to a factor of 1/2, which is due to the fact that this is the contribution of right-movers only, the obtained result is 
just the DOS of free massive particles with dispersion e{q) — \/ v^q^ + A^. Since there ai:e,no particles above energy 



A at zero temperature, there are no interaction corrections to the density of stat£ 
this result consists in using the Ising-model representation of B fields (see Ref.E 



1 as wellEil. Another way of deriving 
for details) . In this representation, 



where [ii (jJi) are order (disorder) fields of the d 
when At ^ 1, 



exp(— i-\/7r0) = (71/12 — i/ii(72, 

exp(j-\/7r(p) ~ fj,ifi2 + i<7icr2, (119) 
2 Ising model, whose correlators are known. At long times, i.e., 



(/i,(T)/X,(0)) 
((7,(t)(7,(0)) 



S,,KoiA\T\). 



(120) 



As a resu lt, Fp{T) ~ sgn{T)KQ{A\T\). Because of the condition At ^ 1, there is no discrepancy between Eqs. ( |l20| ) 
and ( 117 ), since the leading asymptotic term of Kii{x) is j/- independent. Hence one again finds a square-root singularity 
in Fp'^^{ijj) for — A <C A. The correspondence with the Ising model allows one to estimate neglected terms as e~^^'^ . 

A more illuminating way to understand the square-root singularity is provided by the real-time calculation. Starting 
from Eq.(117), it can be shown thatEJ 



ImlFj;"' 



A 

47rw 



/ dte"^' (Ki i-iAt) - Ki (iAt) 

J — 00 ^ 



' —00 



A f°° 

— / dts\\\(ujt)Yi{At) = e(cj - A) 
2w Jo 



A2' 



(121) 



where Yi{x) is the Bessel function of the second kind. Since asymptotically Yi{x) ~ s in(x)/-v/x, the origin of the 
singularity at w = A can be easily understood. For — A <C A, the integrand of ( |121| ) oscillates very slowly, with 
period = 2it/{ijj — A). The integral is thus determined by long times, t to, and can be estimated as Jp" ^ ^/to. 
We see that the threshold behavior of the DOS is determined by times much longer than 1/A, which justifies our use 
of the long-time asymptotics of Bessel functions to evaluate the DOS at w « A. 



3. DOS of a physical electron 



To find the density of states of a physical electron, we have to convolute Eq.(118) with the contribution of the 
gapless spin mode: 



piu;,x^ 0) 



de„ 
2tt 



■Fp{en,X = 0)F„{LOn 



\uJn — — tUJ 



- r deIm[F;^\e)]Im[F:^'{LU-e) 

TI" Jo 



(122) 
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Since Fo-(e„) y^a/iven, we find 



p{uj,x = 0) 



'\l V J A 



1 



For — A ^ A, where our derivation is valid, 



= 0) 



71" /ff 



eiuj-A). 



(123) 



(124) 



Instead of a square-root singularity of Eq. (118), the DOS of a physical electron exhibits a regular behavior approaching 
a finite value at the threshold. This modification is due to dressing of the gapped charge mode by gapless spin 
excitations. At energies much above the gap DOS increases, p{uj,x = 0) ~ ^/uj, which means that the speettal 
weight is shifted to higher energi es. The energy-independent electron DOS near the threshold was obtained inc^"^^ 



Parenthetically, functional form ( 124 ) rem ains valid when the spin channel is gapped as well. In this case the density 
of states of spin excitations is given by Eq.( |ll^ ) with A A^-. Peforming integration in (122), we again find behavior 
described by Eq. (124) near the threshold A + A^. In other words, the gap in the electron's DOS is given by the sum 
of charge- and spin-sector gaps. 



C. DOS of the electron ladder. 



We consider now the tunneling density of states of a two-channel wire. To apply Eq.(lOS), one needs to have 
the strong-coupling fixed-point values of K^,, which are not known for a general case of two non-equivalent channels 
coupled by the Coulomb interaction. Just to illustrate what kind of behavior one might expect in this case, we consider 
an electron ladder with Kii, = K21, = p, cr) in the Cooper phase. The right-moving fermion is represented by 



Rn=l,s 



Now we apply Eq.(108) to the correlator F{t) = —{TrRi^s{T,Q)R[ ^{<d,Q)) . In the Cooper phase, 
gapped, hence their conjugates are exponentially suppressed. As a consequence, e.g.. 



On the other hand, Lpp-^- and 9p+ remain critical. As a result 



F{T)(xsgn{T) \t\ exp - — {2m„Va + 'mp^Vp^ 



\ 1/8 




j exp 


16 



K= \{i + Kp 



yKp 



Hence the DOS behaves as 



ejuj - Asc) 

py^) « 7- — X — v7' 7 



(^ 



-i5~Kp+~l/Kp+), 



(125) 
and Lpa-± are 

(126) 

(127) 
(128) 



where Asc = (7r/16)(2TO„u„ + nip^Vp-). Note that 7 < 3/8 < 1/2 for Kp+ < 1. 

Exactly at half-filling, when p+ mode is also gapped due to Umklapp scattering and there are no more gapless 
modes, our procedure gives F{t) ^ e~'^l'^l/\/M' ™ agreement with recent exact resultl£3. The corresponding DOS is 
that of a free mass ive p article, Phf j^) ^ Q{lo ~ A)/^Jijj — A. 

Comparison of (128) with (|124[) shows that softening of the square- root singularity is pronounced weaker for a 



ladder than than for a Hubbard chain, because a ladder has three gapped and only one gapless mode. As repulsion 
in the p-\- channel becomes stronger, i.e., as -ft'p+ decreases, the singularity becomes weaker and disappears at -ftrp+ — 
(5 — /2 « 0.2. For even smaller Kpj^^ we have p{ijj = Asc) = 0. This behavior though is not very realistic as it 
requires very strong repulsion. For weak repulsion, i.e., when Kp^ « 1, 



p{uj) cx <d{uj - Asc)/{^ 



Asc)'/', 



(129) 



and t he th reshold singularity is still present albeit softened compared to the free-massive-particle case. We note tbsit 
p{ll)) (|128D is similar to the DOS of high unoccupied subbands of a quantum wire, considered recently by BalentsEd. 
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One would expect that the long-range order of the Cooper phase affect tunneling. Indeed one finds that pair 
correlations are determined by the Luttinger parameter of the total charge fluctuations onlyll3, 

(i?i,,(T)Li,_,(T) l\ _MR\.s{'^)) - (e-*^''^+(")e'^''^+(°)) - T-^/i^K,+ )^ (^3Q) 

whereas all other two-particle combinations decay exponentially. Thus, although the single-particle density of states 
is strictly zero, the two-particle one is not. In principle, this effect can be checked experimentally by tunneling into a 
two-channel wire from the superconducting tip - one should observe then a nonzero tunneling current of Cooper pairs. 
Its magnitude, however, will be much smaller than the current in a system of a normal tip and gapless wire, because 
the probability of two-particle (Copper pair) tunneling is exponentially smaller than that of single-particle 

tunnehng (|Tp). 



VII. TUNNELING INTO THE END OF A GAPPED WIRE 



Tunneling into the end of a Luttinger liquid is different from tunneling into the bulk. The reason for this difference 
is the open boundary condition — for the electron wavefunction. For boson modes describing charge and spin 
displacements, this condition means pinning at the baundary. The difference between the edge and bulk tunneling 
was considered first theoretically by Kane and FisheiQ, and has recently been observed in experiments on tuneling 
into carbon nanotubesB. A rigorous treatment of a Luttinger liquii-jwith open boundary conditions, which involves 
re-formulation of the bosonization procedure, can be found in Refs.Eil~E3. 

Suppose now that a two-subband wire is drive n into a CDW state by direct backscattering processes accompanied by 



density adjustment, as described in Sec. Ill A 1 . The relative mode of charge excitations is described by Hamiltonian 



(Eq), in which we put Skp = 0. From the equivalence of Eqs. (101) and (HS) (with Skp — 0), we expect the fixed 



point-value of A'_ in the CDW-phase to be the same as for a half-filled Hubbard chain, i.e., Kl_ = 1/2. The total 



charg e mode {ip+) remains gapless and plays the same role as the spin mode of the Hubbard chain [see Eqs. (124) and 



(128)]: it softens the threshold singularity of the DOS. For lu < Aqdw^ the density of states is equal to zero. Thus if 
the tunneling contact probes the interior of the wire, a gapped behavior is observed. 

However, the DOS at the end of a wire exhibits gapless behavior, i.e, Pend oc Iwl""*, as we will demonstrate in the 
rest of this Section. Consequently, I{V) oc |y|"<i+i for tunneling from a macroscopic (Fermi-liquid contact) into the 
end of the wire, and I(V) cx for tunneling through a barrier located somewhere within the wire. 

This very different behavior of the DOS in the bulk and at the end of the wire can be understood physically for 
tunneling through a barrier located within the wire (cf. Fig. ||). Without the barrier, the CDW is free to slide and 
the conductance is the same as in the absence of any interactions. Squeezing one electron into the middle of the 
wire leads to creation of a soliton-like compression in one of the two modes, and to accompanying it "stretch" in the 
other mode, which requires an energy of the order of the charge gap Acdw- However, such an excitation needs not 
be created when the barrier distorts the uniform profile of the CDW. Indeed, the boundary condition imposed by 
the barrier pins the mode (p- at x = to a value which is different from the one in the bulk, ^-{x oo) — ^/tt/S 
(the latter follows from the minimization of the CDW energy, as illustrated in Fig.||). Therefore, the regular order of 
CDW is already frustrated near the barrier: there is a built-in compression in one of the modes, and depression in 
the other. The electron that tunnels through the barrier arrives into the "stretched" mode. Upon the proper shift of 
both modes, the system arrives into a state with the same energy but with a switched "polarity" of frustration. This 
consideration is true if the barrier is strong enough to destroy the CDW order in its vicinity, i.e., the barrier height 
is larger than Acdw- 



FIG. 6. Illustration of the CDW in the vicinity of a barrier : while the bulk value for ip- is ^7r/8, the barrier forces ip- to 
vanish at a:: = 0. 

We now give a derivation for the current through the barrier. To model the boundary conditions corresponding 
to a wire cut into two semi- infinite pieces, we choose the potential barrier in the form u'(a:;,r_L) = TraWd{x), so that 
Wn=i,2 = TTaW and Winter — [cf. Eqs. (80 H^)]. To find the current in, e.g., the 1st subband, through the barrier 

at a; = 

I (X \im{dtipi{x ^0,t))/e, (131) 

t — *oo 



26 



we need to calculate the rate at which (pi{x = 0, t) increases in the stationary limit in the presence of a potential drop 
~ —eV J dx (9a;(/7+/2Y^)sgn(a;), proportional to the applied volt age V (due to their equivalence, it does not matter 
which of the two modes we use for measuring the current). Eq. ( |l3l| ) reduces the transport problem to that of the 
dynamics of a quantum particle q[t) = ipi{x — 0,i) subject to "damping" by all of the remaining bulk degrees-nf 
freedom, including those of the second mode. Therefore, we can employ methods of dissipative quantum mechanicS'Q 
to solve this problem. 

The effective action of the boundary mode Sef / is obtained by tracing over these remaining degrees of freedom 

e-^"«['l = J V[^i]V[^2] 6[qiT) -^x^ 0,t)] S[Mx = 0,t)] e-^, (132) 



where 

1 



S = dx dr 



+ (133) 



2VF 



(drip-) 



VF 

2Kl 



4/6. 2 
2"'^- 



where =b combinations are defined in ([43|). Performing the integration, we get 

Ses[q] = \[ dr I dr' q{T)JC{r - T')q{T'), 



(134) 



(135) 



which accounts for the influence of the bulk modes exactly. In ( |132| ) we assume a large barrier so that incoherent single 
electron passages through subband 1 determine the current, while subband 2 is fixed by the barrier, cf. (132). Simul- 
taneous contributions from subband 2 represent coherence effects, and are of higher order in the barrier transmission 
coefficient. The dynamics of q is governed by 



So = J dr (wcos2y^q - -^<l) + -S'off, 



(136) 



which includes the applied voltage. Eq. (|136| ) describes a damped quantum particle on a tilted washboard potential. 
Tunneling between a djace nt minima of the potential corresponds to single electron transfers. 

For large Eq. (136) maps onto a tight binding model with nearest neighbor hopping T, whose value can in 
principle be deduced from W , assuming a (5-barrier : it is renormalized lw-|damping compared to the bare valueuS. In 
this analogy, (131) corresponds to the particle's mobility, studied in Ref.E3. In leading order ^ |Tp, the result is 



where 



I{V) = e\T\ 



wit) 



dt sm(eVt) Im e 



-w{t) 



J{UJ) 



(1 - e-'"^^) 



10 ^ 

The spectral function J{u}) is related to the Fourier transform /C(cli„) of the kernelEII appearing in (|13 



J(^) 



lim Im /C( 

77^0 



-ri) 



with uJq 



fbsVF /tto^- This yields 



liV) ^ V 



k1 



e(w - Wo), 



(137) 



(138) 



(139) 



(140) 



which gives ad = 1/2K+ — 1 at voltages V < uJo/e. Above the gap, i.e, for V 3> ojo/e, gapless behavior of a 
two-subband Luttinger liquid is restored I{V ^ wg/e) ^ . 

The gapless behavior of fermion's DOS at the end of the wire can be interpreted in terms of a midgap state in 
the (— ) channel localized near the boundary. This is a bound state with zero energy formed in the potential well 
created by the static distortion of the (p- field subject to an open boundary condition (x = 0) = (see Appendix 



27 



D of the paper by Fabrizio and Gogolina). Qualitatively, this state is similar to the boundary state in a doped 
two- leg spin ladder, which represents free S — \/2 spin induced by a charged impurity, see, e.g.,L3. As a result, local 
density of states of the (— ) channel at the end of the wire takes the form p-{e) ^ A5(e) + p!!;^, where p"^* stands 
for contribution of massive modes with energies above the CDW gap. On the other hand, the (+) mode remains 

gapless, and its end-chain DOS is given by p+{x = 0,a;) ~ ijJ^'^+ ^, BEiS. The factor of 1/2 in the exponent is due 
to "factorization" of the electron operator into (±) modes. Another consequence of s uch factorization is that the DOS 
of a physical electron is a convolution of the DOS of the (±) channels [cf. Eq. (|122|)] 



Pphys{x ^0,uj) I dep+(a; = 0,e)p_(a: = 0,w - e). (141) 
Jo 

Therefore, 

Pphysix = 0, w) - p+{x = 0, CJ) = (142) 

which implies = 1/{2K+) — 1, in agreement with the result of the explicit calculation of I{V) presented above. 



Eq. (140) 



VIII. EXPERIMENTAL CONSEQUENCIES AND CONCLUSIONS 



Apart from the rather well-known Luttingcr-liquid phase, a two-subband quantum wire may also exhibit either 
a CDW or superconducting (Cooper) phase. In both of these phases, certain modes of inter-subband charge- and 
spin-excitations are gapped, whereas the center-of-mass charge mode remains gapless. As a result, the conductance 
remains at the universal value of 2e^ /h per occupied subband, irrespective of whether the wire is in a gapless or 
gapped phase. However, the single-particle density of states in the middle of the wire has a hard gap. Above the gap, 
the DOS exhibits a non-universal threshold behavior p{ui) ~ 8(u; — A)(ti; — A)"'', where h < 1/2. Softening of the 
threshold singularity is due to "dressing" of gapped modes by the remaining gapless one. 

We find that the DOS for tunneling into the end of a wire in the CDW phase remains gapless, with the exponent 
determined by the center-of-mass mode only. This effect is due to frustration introduced into the CDW order by an 
open boundary (strong barrier). _ . 

Where should one look for such exotic phases of a quantum wire? We believe that quantum wiresEll prepared by 
cleaved edge overgrowth technique may be well-suited for observing the CDW phase. Indeed, the cross-section of such 
a wire is close to a square, which implies that the lowest states of transverse quantization should be close in energy. 
Hence, one-dimensional subbands can have close Fermi-momenta, which is a necessary condition for the formation of 
the CDW state. The Cooper phase, on the other hand, requires effective attraction in the relative charge channel, 
and has the best chance to occur when the second (upper) subband just opens for conduction, i.e., is near the Van 
Hove singularity. 
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